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We discuss the role of tangential stabilization in a curvature driven flow of planar curves. The governing system of nonlinear 
parabolic equations includes a nontrivial tangential velocity functional yielding a uniform redistribution of grid points along 
the evolving family of curves preventing numerically computed curves from forming various instabilities. 
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CN ■ 1 Introduction 
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O . In this paper we study evolution of a family of closed smooth plane curves Tt : S — > K , t > 0, driven by the normal velocity 
■ V which is assumed to be a function of the curvature k, tangential angle u and position vector x €Tt, 

OO ■ V = I3(x,k,iy) . (1) 

Csl ■ y ^ / 

, As a typical example one can consider a normal velocity of the form; v ~ k (mean curvature driven flow), v = k^ (affine 

invariant flow), v = a{x,v)k + c{x, v) (Gibbs-Thomson law), etc. Geometric equations of the form ^ can often be found 
in variety of applied problems like e.g. the material science, dynamics of phase boundaries in thermomechanics, in modeling 
of flame front propagation, in combustion, in computations of first arrival times of seismic waves, in computational geometry, 
f~| , robotics, semiconductors industry, etc. They also have a special conceptual importance in image processing and computer 
vision. For an overview of important applications of ([T]i we refer to a book by Sethian [1]. 

An idea behind the direct (or Langrangean) approach consists in representing the family of immersed curves Tt by the 
position vector x S M^, i.e. Vt = Image(a;(., t)) = {x{u, t), u € S^} where a; is a solution to the geometric equation 



< 
7^ 



^ . dtx = 13N + aT (2) 

> ' 

■ where (3 = I3{x,k,v), N ~ (— sin i^, cos and T ~ (cos i/, sin are the unit inward normal and tangent vectors, respec- 
T— I , tively. We chose the orientation of the tangent vector T such that det(T, iV) = 1. Notice that the presence of arbitrary 
CO ■ tangential velocity functional a has no impact on the shape of evolving curves and thus a can be viewed as free parameter to 

, be suitably determined. The unit arc-length parameterization of a curve T = Image(x) will be denoted by s. Then ds ~ g du 

' where g = \dux\. 

According to [2, 3] (see also [4, 5]) the system of governing equations for the curvature fc, tangent angle local length g 
and the position vector x reads as follows: 



o 



>; dtk^d^P+adsk+k^P, dtv ^ (3'kdlv^{a+[3l)dsV+\/ ^p.f , dtg ^ -gkf3+dua , dtx ^ (3N+af (3) 

where (u, t) ^ x (0, T), ds = gdu, A solution (fc, v, g, x) to ^ is subject to initial conditions and periodic boundary 
5—1 conditions in the u variable. 

a ' 

_ _ I 

2 The role of the tangential velocity functional 

Notice that the functional a is still undetermined and it may depend on variables k, i', g, x in various ways including nonlocal 
dependence in particular. Suitable choices of the tangential velocity functional a are discussed in a more detail in this section. 
Although a plays an important role in the governing equations resulting in dependence of k, g, x on a, the family of planar 
curves Tt = Image(a;(., <)), t e [0, T), is independent of a particular choice of a. 

To motivate further discussion, we recall some of computational examples in which the usual choice a = Q fails and may 
lead to serious numerical instabilities like e.g. formation of so-called swallow tails. In Figure [T]-a) we computed the mean 
curvature flow of an initial curve (bold faced curve). We chose a = 0. It should be obvious that numerically computed grid 
points merge in some parts of the curve Tt preventing thus numerical approximation of Ft , t G [0 , T) , to be continued beyond 
some time T which is still far away from the maximal time of existence Tmax- This and many other examples from [2, 3] 
showed that a suitable grid points redistribution governed by a nontrivial tangential velocity functional a is needed in order to 
compute the solution over its life-span. 
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Fig. 1 Numerically computed evolution of an initial curve (bold); a) merging of numerically computed grid points in the case of zero 
tangential velocity q = 0; b) enhancement of grid point redistribution in the case of asymptotically uniform parametrization. 



The idea behind construction of a suitable tangential velocity functional a is rather simple and consists in the analysis of 
the quantity 9 defined as 6' = ln{g/L) where g ~ \dux\ is a local length and L ~ Lt = /p ds = g{u, t) du is a total length 
of the curve Fj = Image(a;(., <)). The quantity 6 can be viewed as the logarithm of the relative local length ratio g/L. Taking 
into account equations Q and the equation for the total length j^L + J^, k(3ds = (obtained again from ^ by integration) 
we have 

dt0 + kf3 - {kp)r ^ dsa (4) 

where (fc/3)r denotes the average of over the curve F, i.e. (fc/3)r = x /r ^P*^^ an appropriate choice of dga in the 
right hand side of ^ appropriately we can therefore control the behavior of 6. Equation (|4| can be also viewed as a kind of a 
constitutive relation determining redistribution of grid point along a curve. The simplest possible choice of dsa is: 

dsa ^kf3- {k(3)r (5) 

yielding dtO = in (UJi. Consequently, g{u,t)/Lt ~ g{u,0)/Lo for any u ^ S^, t ^ [0,Tmax)- Notice that a can be 
uniquely computed from ^ under the additional renormalization constraint: a{0, t) = 0. The tangential redistribution driven 
by a solution a to (|5]l is refereed to as a parameterization preserving relative local length (c.f. [2]). It has been first discovered 
and utilized by Hou et al. in [6, 7] and independently by the authors in [2, 3]. 
A more general choice of a is based on the following setup: 

d,a = kf3- {k(3)r + (e"^ - l) Lu{t) (6) 

where u) G Ll^^{[0, Tmax))- If we additionally suppose /q^™""^ ^{t) = +oo then, after insertion of Q into ^ and solving 
the ODE dtO = (e^^ — l) Lo{t), we obtain 9{u,t) ^ as t Tmax and hence g{u,t)/Lt — > 1 as t ^ T,nax uniformly 
w.r. to u e S^. In this case redistribution of grid points along a curve becomes uniform as t approaches the maximal time of 
existence Tmax- We will refer to the parameterization based on ^ to as an asymptotically uniform parameterization (c.f. [3]). 
The impact of a tangential velocity functional defined as in (|5]i on enhancement of redistribution of grid points can be observed 
from two examples shown in Fig. [l]-b) computed by the authors in [2]. It can be shown that the appropriate choice for the 
control function uj takes the form uj = ki + H2{k(3)r and ki,K2 > are given constants. A detailed discussion on this 
topic can be found in [3,4]. If we insert tangential velocity functional a computed from (|6]l into (O the system of governing 
equations can be rewritten as follows: 

dtk = dlP + ds{ak) + fc(fc/3)r + (1 - L/g) kuj , dtiy = Pkd^v + (a + /3^)5,i^ + V^/3.f, (7) 

dtg = -9{k(3)r + (i - g)uj , dtx ^ /3N + af . 
It is worth to note that the strong reaction term in (|3]l has been replaced by the averaged term fc(fc/3)r in (|7]i. This is a very 
important feature as it allows for construction of an efficient and stable numerical scheme discussed in more details in [3-5]. 
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